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Critical behavior of the Ising model is investigated at the center of large scale finite size systems, 
where the lattice is represented as the tiling of pentagons. The system is on the hyperbolic plane, 
and the recursive structure of the lattice makes it possible to apply the corner transfer matrix 
renormalization group method. From the calculated nearest neighbor spin correlation function 
and the spontaneous magnetization, it is concluded that the phase transition of this model is 
mean-field like. One parameter deformation of the corner Hamiltonian on the hyperbolic plane 
is discussed. 
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§1. Introduction 

Baxter's method of corner transfer matrix (CTM) has 
been known as one of the representative tool for analyti- 
cal stydy of statistical models in two dimension (2D). 1-3 ) 
The method is also of use for numerical calculations of 
one point functions, such as the local energy and the 
magnetization. 3 ' This numerical application is a kind of 
numerical renormalization group (RG) method, where 
the block spin transformation is obtained from the di- 
agonalization of CTMs. Such a RG scheme has many 
aspects in common with the density matrix renormal- 
ization group (DMRG) method, 4 ~ 7 ' expecially when the 
method is applied to 2D classical lattice models. 8 ' 

Introducing the flexibility in the system extension pro- 
cess of the DMRG method to the Baxter's method of 
CTM, the authors developed the corner transfer matrix 
renormalization group (CTMRG) method. 9 ~ 12 ' In this 
article we report a modification of the CTMRG method, 
for the purpose of applying the method to classical lat- 
tice models on the hyperbolic plane. Using the recursive 
structure of the lattice, we obtain one point functions at 
the center of sufficiently large finite size systems. 

Quite recently Hasegawa, Sakaniwa, and Shima re- 
ported deviations of critical indices of the Ising model on 
the hyperbolic plane from the well known Ising univer- 
sality classes in two dimension. 14, 15 ' They predicted that 
phase transition of such systems would be mean-field like. 
To confirm their prediction, in the next section we con- 
sider the Ising model on a lattice, which is represented 
as the tiling of pentagons. The necessary modification of 
the CTMRG method on this lattice is explained in §3. 
We calculate the nearest neighbor spin correlation func- 
tion and the spontaneous magnetization at the center of 
large scale finite size systems. Critical indices for these 
one point functions are studied in §4, and we confirme the 
mean-field like properties of the phase transition. Con- 
clusions are summarized in the last section. We discuss 
a possible deformation of corner Hamiltonian in the hy- 



perbolic plane. 

§2. Ising Model on the Tiling of Pentagons 

Let us consider the hyperbolic plane, which is the two 
dimensional surface with constant negative curvature. 
Figure 1 shows a part of a sufficiently large regular lat- 
tice on the plane, where the lattice is constructed as the 
tiling of pentagons. 13 ' All the arcs are geodesies, which 
divides the lattice into two parts of similar structure. 
Each lattice point represented by an open circle is the 
crossing points of two geodesies. 




Fig. 1. Ising model on a regular lattice on the hyperbolic plane. 
Open circles represent Ising spins on the lattice point, and 
W(a a , o~, , <r c , <r ,, <r e ) is the local Boltzmann weight denned in 
Eq. (2.2). Two geodesies drawn by thick arcs divide the system 
into four quadrants, which are called as the corner. 



We investigate the ferromagnetic Ising model on this 
lattice. The Hamiltonian of the system is defined as the 
sum of nearest neighbor Ising interactions 



H 



(2.1) 



where represents pair of neighboring sites, and 

a i — ±1 and a- = ±1 are the Ising spins on the lattice 
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points. Throughout this article we assume the absence 
of external magnetic field. For the latter conveniences, 
we represent the system as the 'interaction round a face 
(IRF)' model. The local Boltzmann weight for each face 
of pentagonal shape — the IRF weight — is given by 



W(a a ,a b ,a c ,a d ,a, 
(3J 



(2.2) 



= exp 



a a a b + a b a c + a c a d + a d°e + °e°a, 



where (3 = l/k B T is the inverse temperature, and a a , a b , 
<r c , <r d , and a e are the spin variables around the face as 
shwon in Fig. 1. 

The partition function of a finite size system (with suf- 
ficiently large diameter) is formally written as the con- 
figuration sum of the Boltzmann weight of the whole 
system 



Z =J2 U w 



(2.3) 



all the all the 
spills faces 



where we are interested in the thermodynamic limit of 
this system. Note that it is rather hard to investigate the 
system by use of the Monte Carlo simulations, since the 
number of sites contained in a cluster blows up exponen- 
tially with respect to its diameter. As a complcmental 
numerical tool, we employ the CTMRG method. 

§3. Corner Transfer Matrix Renormalization 
Group Method 

The two geodesies shown by thick arcs in Fig. 1 divide 
the system into four parts, which are called as corners. ^ 
Figure 2 shows the structure of a corner. We label the 
spins on a cut as {a 1 , a 2 , <t 3 , . . .}, and those on another 
cut as {<j[, <7 2 , erg, . . .}, where <r 1 is equivalent to (t[. The 
corner transfer matrix is the Boltzmann weight with re- 
spect to a corner, which is calculated as a partial sum of 
the product of IRF weights in the corner 



C(<?i, 02, 03, • • • I ^i,cr 2 ,cr 3 , . . .) 

= e n w - 



(3.1) 



spins inside faces in 

the quadrant the quadrant 



The configuration sum is taken over spins 'inside' 
the corner, leaving those spins on the cuts. Con- 
ventionally the matrix C is interpreted as block di- 
agonal with respect to cr l and a[, and the element 
C(a' 1 , 2 , 3 , • • • | (T 1 , ct 2 , 03) ■ • ■) f° r those cases a 1 ^ cr[ 
is set to zero. The CTM thus defined is symmetric in 
the case of the pentagonal lattice under consideration. 

As shown in Fig. 2 a corner has the structure where 
three parts labeled by P and two parts labeled by C 
are joined to a face W . The 'fusion' relation can be 
represented by a formal equation 9, 10 ) 



C = W ■ PC PC PC . 



(3.2) 



Note that C is a corner of smaller size. 

For convenience, let us observe the structure of the 
part of the system shown in Fig. 3, where two P, and C 
are joined to W. Labeling the shown part by P, we can 




Fig. 2. Recursive structure of a corner C. 

formally write the fusion relation in the same manner 

P = WPCP. (3.3) 

For a conventional reason we call P as the 'half-row', 
although P is not a row on the hyperbolic plane. It is 
easily understood that P is a half-row of smaller size. 
We have thus obtained recursive structure of the corner 
C and the half-row P. 




Fig. 3. Recursive structure of the 'half-raw' P. 



As we have defined CTM for a corner, let us express 
the Boltzmann weight with respect to the half-row 



P(a 1 ,a 2 ,a 3 ,...\a 1 ,a 2 ,a 3 ,...) 

= e n w 



(3.4) 



the half-row the half-] 



in the matrix form, where the positions of spins 
{0 n , 2 , ct 3 , . . .} and {a[, a' 2 , cr' 3 , . . .} are shown in Fig. 3. 
We call the weight in the matrix form as the half-row 
transfer matrix (HRTM). 
The 4-th power of the CTM 



p = C A 



(3.5) 



is a kind of density matrix, since its trace gives the par- 
tition function 

Z = T r( o = TrC 4 (3.6) 
of a finite size cluster that consists of four corners. The 
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matrix dimension of CTM , and also that of the HRTM, 
increases exponentially with respect to the system size. 
In order to obtain Z numerically up to sufficiently large 
systems, we introduce the block spin transformation that 
is created from the diagonalization of the density matrix 

p.4,5) 

Instead of directly diagonalizing p in Eq. (3.5), we first 
create its contraction 

P'(°2>°3 ■■■\ a 2>< T 3i-- •) 

= X! P( CT 1> £7 2> CT 3---| CT 1>°'2> <7 3>---) ( 3 - 7 ) 
a' 1 =<7 1 =±l 

and then diagonalize it 

p'(cr' 2 ,a' 3 ... | <7 2 ,0- 3 , . . .) 

= ^^,4...|0A^(<7 2 ,<7 3 ,...|0, (3.8) 

where the eigenvalue A^ is non-negative. Following the 
convention in DMRG, we assume the decreasing order 
for A^. The orthogonal matrix A(a 2 , tr 3 , . . . | £) repre- 
sents the block spin transformation from the 'row-spin' 
{cr 2 , (T 3 , . . .} to the effective spin variable £. We keep m 
numbers of representative states, which correspond to 
major eigenvalues, for the block spin variable £. Apply- 
ing the matrix A to the CTM and the HRTM, we obtain 
'renormalized matrices' of 2m-dimension 

C(ct[, a 2 , ...\a 1 ,a 2 ,...)^ C(a',£'\a,£) 

P(a[,a 2 , . . . | a,, a 2 , . . .) - P(a', ? \a,£), (3.9) 

where we have dropped the indices from a 1 and er^. 9,10 ) 
Combining the recursive structures in Eqs. (3.2) and 
(3.3), and the renormalization scheme in Eqs. (3.7)-(3.9), 
we can obtain the CTM and HRTM in the renormalized 
form for arbitrary system size by way of successive ex- 
tension of the system. 




(a) (b) 

Fig. 4. Extension pocess of (a) CTM and (b) HRTM in the renor- 
malized expression. 



Suppose that we have C(cr', £' | er, £) and P(cr', £' | er, £) 
for a finite size cluster. In order to explain the extension 
process, let us rewrite these matrices as C(s' , (' \ s, and 
P(s',C\s,(). 

(1) Substitute C(s',C\s,0 and P(s',C\s,0 into the 
fusion process in Eqs. (3.2) and (3.3). Figure 4 
shows these fusion processes among W, C, and P, 
where rectangles correspond to the block spin vari- 



ables. The spin variables that are contracted out are 
shown by black marks. As a result, we obtain the 
extended CTM C(er' , s' , (' \ a, s, C) and the extended 
HRTM P(a' ,s' ,C\v,s,()- 

(2) From the extended CTM C(er' , s' , (' \ a, s, C) obtain 
the density matrix p(a' , s' , ('\ <r, s, () by Eq. (3.5). 
Contracting out the spin at the center, obtain 
p' (s' , £' | s, C) as Eq. (3.7), and diagonalizing it to ob- 
tain the block spin transformation matrix A(s, £ ^) 
from Eq. (3.8). 

(3) Applying A(s, C| to both C(a' , s' , C'| o, s, () and 
P(a',s',C\ c,s,C), obtain the extended CTAI 
C (a 1 , £' | er, £) in the initial form, and the same for 
HRTM to obtain P(a' , £'| a, £). 

(4) return to the first step. 

The system size, which is the length of the longest 
geodesies in the system, increases by 2 for each itera- 
tion. 16 ) In order to start the above extension process, we 
set the initial condition 

C(a' | a) = P(a' \ a) = 5(a' \ 1) 6(a\ 1) (3.10) 

that represents ferromagnetic boundary, where 5(a\b) = 
S a h is the Cronecker's delta. 

During the iteration we can obtain one point functions 
at the center of the system. For example, the sponta- 
neous magnetization is calculated as 

X a P( a > s ' CI CT ' s ' o 



Tp p( a > s ' ci a > s > o 

In the same manner we obtain the nearest neighbor spin 
correlation function (<rs). It should be noted that one 
point functions thus calculated at the center do not al- 
ways represent the averaged property of the whole sys- 
tem even in the thermodynamic limit, since the area near 
the boundary has non-negligible weight in the hyperbolic 
plane. 

§4. Numerical Result Compared With the 
Bethe Approximation 

Let us calculate the spontaneous magnetization (er), 
and spin correlation function (as) for the nearest spin 
pair. We regard the Ising interaction strength J as the 
energy unit, and use the temperature where the Boltz- 
mann constant fc B is equal to unity. Most of the numer- 
ical calculations are performed keeping m = 40 states. 
The dumping of the density matrix eigenvalues is very 
fast, and actually the calculated results with m = 10 do 
not differ from those obtained with m = 40 even at the 
critical temperature T c . The iteration number required 
for the numerical convergence is at most 400000 for the 
calculated data points. 

Figure 5 shows the calculated results. The square of 
the spontaneous magnetization (a) 2 is a linear function 
of temperature in the neighborhood of T c . From the 
behavior we estimate the transition temperature T c — 
2.799. The nearest neighbor spin correlation function 
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Fig. 5. Square of the spontaneous magnetization (a) 2 (upper) 
and the nearest neighbor spin correlation function (as) (lower) 
with respect to the temperature T. 



(as) has a kink at T c , and is linear in T around there. 
These calculated results support the existence of mean- 
field like transition, that is subject to the critical indices 
(3=1/2 and a = 0. We thus confirmed the prediction 
by Hasegawa, Sakaniwa, and Shima. 14 ' 15 ' 

Compared with the transition temperature of the 
square lattice Ising model ySq uarc = 2.269, the calcu- 
lated T c is fairly higher and is close to the transition 
temperature calculated from the Bethe approximation 
T Bctho = 2. 885 . 17,18) The result SU gge S t that neglection 

of the 'loop back effect' is not so conspicuous in the hy- 
perbolic plane. 

§5. Conclusion and discussion 

We have calculated the spontaneous magnetization 
and the nearest neighbor spin correlation function of the 
Ising model on a pentagonal lattice on the hyperbolic 
plane. The numerical algorithm of the CTMRG method 
is modified for this purpose. The calculated critical tem- 
perature is T c = 2.799, and we observe the mean-field 
like phase transition. 

The modified CTMRG method we have developped is 
applicable to regular lattices that consists of geodesies 
on the hyperbolic plane. For those lattices that does 
not contain geodesies, one has to either treating asym- 
metric density matrix or to draw geodesies by use of 
transformation such as duality transformation and the 
star-triangle relation. Generalization of the modified 
CTMRG method to the vertex model is straight for- 
ward. It may be interesting to classify ordered states 
of eight-vertex model on a variety of regular lattices in 
the hyperbolic plane. 

An interest is in the eigenvalue structure of the den- 
sity matrix at the transition temperature. Its analytic 
form is not well defined in the thermodynamic limit of 
classical lattice models on the flat 2D plane. 19,20 - 1 The 
rapid eigenvalue dumping observed on the hyperbolic 
plane suggests that there would be a way of regulariz- 
ing the CTM at the criticality. The classical-quantum 



correspondence from such a view point is worth consid- 
ering. Formally speaking the corner transfer matrix C 
can be written as the exponential of the corner Hamil- 
tonian H c . When the lattice is on the flat plane, the 
simplest example of the corner Hamiltonian is written in 
the sum of local operators 



H c = h{a 1 ,a 2 ) + 2h(a 2l cr 3 ) + 3h{a 3 ,a i ) 
= h 1 + 2h 2 + 3h 3 + 4h 4 + ... , 



(5.1) 



where h i = h(a i ,ij i+1 ) is the local hamiltonian that acs 
between neighboring sites. A possible one parameter de- 
formation of the corner Hamiltonian to the hyperbolic 
geometry may given by 

/ » \ , sinh2A, sinh3A, 
H C (A) =h,+ . - ; h 2 + . - ; h 3 + . . . , (5.2 
cy ' 1 sinhA 2 sinhA 3 ' v ' 

which is reduced to H c in Eq. (5.1) in the limit A — > 0. 
The deformed Hamiltonian satisfies the recursive struc- 
ture 

sinh 3A , 



H C (A) = cosh A h 



sinh2A^ 
sinh A 3 



sinh A 



+ h 1 + cosh Ah 2 + cosh 2Ah 3 + . . . 



(5.3) 



discussed in Okunishi's RG scheme on the corner Hamil- 
tonian. 21 ) Such a deformation has similar regularization 
effect proposed by Okunishi quite recently 22 ) 
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